Theory of Magneto-resistance of Disordered Superconducting Films 
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Recent experimental studies of magneto-resistance in disordered superconducting thin films reveal a 
huge peak (about 5 orders of magnitude) . While it may be expected that magnetic field destroys su- 
perconductivity, leading to an enhanced resistance, attenuation of the resistance at higher magnetic 
fields is surprising. We propose a model which accounts for the experimental results in the entire 
range of magnetic fields, based on the formation of superconducting islands due to fluctuations in 
the superconducting order parameter amplitude. At strong magnetic fields Coulomb blockade in 
these islands gives rise to negative magneto-resistance. As the magnetic field is reduced the effect 
of Coulomb blockade diminishes and eventually the magneto-resistance changes sign. Numerical 
calculations show good qualitative agreement with experimental data. 

PACS numbers: 71.30.+h,73.43.-f,73.43.Nq,74.20.Mn 



The interplay between superconductivity and disorder 
is a long-standing problem, dating back to the late fifties 
P, 0] , resulting in the common wisdom that weak disor- 
der has no dramatic effect on superconductivity. Strong 
disorder, however, may have a profound effect, driving 
the system from a superconducting (SC) to an insulating 
state. Such a SC-insulator transition (SIT) was observed 
in two-dimensional amorphous superconducting films Q . 
Reducing the film thickness or an increasing perpendicu- 
lar magnetic field drives these films (which are held below 
their bulk critical temperature) from a SC state, charac- 
terized by a vanishing resistance as T — > T c , to an in- 
sulating state, characterized by a diverging resistance as 
T — » 0. The possibility of tuning the system continu- 
ously between these two phases is a manifestation of a 
quantum phase transition 0. 

The origin of this transition is still in debate. While 
some theories |5j claim that it may be understood in 
terms of Cooper-pair scattering out of the SC condensate 
into a Bose-glass state (so-called "dirty boson" models), 
there is evidence, both experimental [a, Li la l2| and the- 
oretical 0| ) that a percolation description of the SIT is 
more adequate for, at least, some of these samples. 

More insight into the nature of the SIT may be gained 
by looking at the magneto-resistance (MR) on its in- 
sulating side. Decade-old experiments observed non- 
monotonic MR, exhibiting a shallow peak at some mag- 
netic field -Bmax 0, E| • Recent experiments ^jj show, 
however, that in some samples the effect is dramatic, 
with the resistance value at the peak, i? ma x7 reaching 
as high as a few orders of magnitude its value at the 
SIT. As the magnetic field is further increased the MR 
drops back a few orders of magnitude (inset of Fig. 3). 
Further investigations of this effect Q reveal that disor- 
der also has a major influence. With increasing disorder 
strength, the critical field for the SIT, B Cl diminishes, 
while i?max increases. The temperature dependence of 
the MR at high temperatures fits an activation-like be- 



havior, R oc exp(To/T), with a magnetic field dependent 
To as seen in the inset of Fig. 4(a). At lower tempera- 
tures there is a deviation from this behavior towards some 
weaker temperature dependence. While an enhancement 
of the resistance with increasing magnetic field may be 
understood in terms of decreasing SC order, the suppres- 
sion of the resistance with further increase of magnetic 
field remains a puzzle. 

Here we propose a model for the MR in the entire range 
of magnetic fields. The model is based on three assump- 
tions. The first is that disorder induces formation of SC 
islands (SCIs) due to fluctuations in the amplitude of the 
SC order parameter. This concept has already been used 
to interpret some experiments [ill ll3T |. It was dire ctly 
observed in other systems by STM measurements Il5 l 
and further corroborated by numerical simulations [lfij . 
The second assumption is that as the magnetic field is 
increased, the concentration and size of these SCIs de- 
crease. Preliminary numerical results support this pic- 
ture and will be reported in the future. The third as- 
sumption is that the SCIs have a charging energy, and 
thus, a Cooper pair entering the SCI via an Andreev tun- 
neling process, has to overcome the charging energy. The 
charging energy is expected to be inversely proportional 
to the island size, and thus to increase with increasing 
magnetic field. 

In order to see the mechanism by which the MR can be 
negative, consider such a system in the strong magnetic 
field regime, B » B max (Fig. 1(a)). Due to the strong 
magnetic field the SCIs are small and have a large charg- 
ing energy. There are two types of trajectories available 
for electron transport: those which follow normal areas of 
the sample ("normal paths", solid lines in Fig. 1(a)) and 
those in which an electron tunnels into a SCI via the An- 
dreev channel ("island paths", dashed lines in Fig. 1(a)). 
The resistance of the normal paths has some value (which 
may depend on e.g. length, temperature, etc.) and is 
assumed to be only weakly affected by magnetic field. 
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Due to the Coulomb blockade, transport through the SCI 
paths is thermally activated, and hence the resistance of 
the island paths is of the form R ~ exp(E c /T), where E c 
is the charging energy of the island. If E c is large then the 
main contribution to the conductance is due to transport 
along the normal paths. Consistent with experiment, the 
MR in this regime is small. 

As the magnetic field is decreased (but still in the 
regime B > B max , Fig. l.(b)), more SCIs are created and 
their size increases, but they are still small enough such 
that transport along normal paths is favorable. How- 
ever, some paths which were normal at higher fields 
(e.g. bottom solid line in Fig. 1(a)) now become island 
paths and hence unavailable for electron transport (bot- 
tom dashed line in Fig. 1(b)). Thus, the effective phase 
space available for electron transport diminishes, result- 
ing in a negative MR. Eventually, at a certain magnetic 
field B — £?max (Fig. 1(c)) some SCIs are large enough 
so that their charging energy is small and the resistance 
through them is comparable to the resistance through 
normal paths. At this point the resistance reaches its 
maximum value, since as the magnetic field is further de- 
creased (Fig. 1(d)) the SCIs are so large that transport 
through them is always preferred over transport through 
normal paths. Increase in number and size of the SCIs 
will thus result in a decrease in the resistance. At the 
critical field B c the SCIs percolate through the system, 
resulting in an insulator-to-SC transition. We note that 
if the temperature is smaller than the Josephson cou- 
pling between two islands then SC correlations will ex- 
tend between them, effectively joining them to a single 
island. Thus, we expect that at low temperatures the 
geometry will be temperature-dependent, resulting in a 
temperature-dependent critical point as is indeed seen in 
some experiments. 

In order to substantiate these heuristic arguments we 
model the thin film by a square lattice, where each site 
can be either normal, with probability p, or SC, with 
probability 1 — p, corresponding to the concentration of 
normal areas and SCIs in the sample. The probability 
p is a function of the magnetic field B, and we may as- 
sume that p(B) is an increasing monotonia function. To 
describe the disordered system in stron g m agnetic fields, 
in the absence of SCIs, we follow Ref. [l?J and assign a 
resistance between any two normal sites of the form 
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where Rq is a constant, is the distance between sites 
i and j, £/ oc is the localization length, is the energy of 
the i-th site measured from the chemical potential (taken 
from a uniform distribution [— W/2, W/2]) and T is the 
temperature. The localization length £/ oc is taken to be 
small (in units of lattice constant), effectively allowing 
only nearest- and next-nearest-neighbor hopping. All SC 
sites that are linked to each other are considered a single 
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FIG. 1: Schematic representation of the model :(a) At strong mag- 
netic fields, B >> -Bmax, The system is composed of small super- 
conducting islands with large charging energy. In this regime trans- 
port through normal paths (solid lines) is always preferable than 
transport through the superconducting islands (SCIs) (see text), 
(b) As the magnetic field is decreased, but still B > B 

max , more 

SCIs appear, resulting in a decrease in available trajectories for 
transport (bottom solid line in (a) and bottom dashed line in (b)) 
and hence negative magneto-resistance, (c) At a certain field -B max 
the resistance of normal paths and paths that include SCIs is com- 
parable, resulting in a peak in the resistance, (d) For even lower 
fields, B c < B < B max, transport through SCIs is always favored, 
resulting in positive magneto-resistance. 



SCI. The resistance between two (neighboring) SC sites, 
i?ssj is taken to be very small compared with Eq.pJ, but 
still not zero and temperature dependent, in such a way 
that it vanishes as T — > (distant SC sites are discon- 
nected). The calculations of resistance were conducted 
with several functional forms for Rss(T) (power law, ex- 
ponential dependence, etc.) and no qualitative difference 
between them was found. The resistance between a nor- 
mal site and a SC site (local N-S junction) is taken to be 



R NS oc exp(£; c /fcT), 



(2) 



where E c is the charging energy of the island. For sim- 
plicity, and to avoid additional parameters, the charging 
energies of all islands were taken to be the same, inde- 
pendent of island size. It is demonstrated below that the 
main experimental observations can be well understood 
even under such an assumption |l8| . 

The lattice is connected to electrodes (Fig. 2) and 
the resistance is calculated numerically using Kirchhoff 's 
laws. The left-most and right-most links in Fig. 2 are 
taken to be SC, thus avoiding a strong dependence of 
the resistance on the properties of the edge sites. In 
Fig. 3 we plot the calculated resistance (on a log scale) 
as a function of the probability p for different temper- 
atures. The calculations were conducted on lattices of 
size 25 x 25 (and repeated for different lattice sizes, with 
no qualitative difference in the results) and log-averaged 
over 100 realizations of disorder. A peak in the MR is 
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FIG. 2: The lattice model is composed of regular sites and super- 
conducting sites. Clusters of the latter form SCIs (shaded islands). 
The resistance between two normal sites (small circles) is given by 
Eq. The resistance between two SC sites (filled squares) is very 
small (and becomes smaller with decreasing temperature, see text) 
for neighboring sites (wavy lines) and is infinite (no link) between 
non-neighboring sites. The resistance between neighboring normal 
and SC sites (thick lines) is much higher than the resistance be- 
tween two normal sites, and is exponential in the charging energy 
of the SC island (Eq.|2J. 

observed at p max = 0.5, with peak resistance four orders 
of magnitude larger than the resistance at the transition. 
The results are compared with the experimental data of 
fl3| (inset of Fig. 3), and good qualitative agreement is 
evident. Notice that the critical probability, p c , defined 
as the probability at which the resistance is temperature- 
independent, is shifted from the percolation critical prob- 
ability. The reason for this is that the resistance of the SC 
links, Rss, is finite. As the temperature is decreased the 
critical probability moves towards the percolation critical 
probability, eventually reaching it at T = 0. 
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FIG. 3: Numerical results : resistance (on a log scale) as a function 
of probability p for different temperatures, with the parameters 
W = 0.4, E c = 4, g ioc = 0.1, T = 0.1, 0.2,..., 4. This is to 
be (qualitatively) compared with the experimental data (inset) of 

m. 

When fitting the resistance as a function of tempera- 
ture to an activation-like behavior, R oc exp(To/T), we 
find a non-monotonic dependence of To on the probabil- 
ity p (Fig. 4), resembling the experimental data of 0] 



(upper inset of Fig. 4) and of [14J. The activation tem- 
perature rises from Tq w W at p — 1 to To w E c for 
p = p max . It then drops back again due to increasing 
weight of SC areas in the sample, eventually reaching 
To = at the transition. Thus, in our picture To is de- 
termined by the Coulomb blockade energy and not by 
the SC gap, as was initially suggested in 13] . 

At high temperatures the activation fit is excellent for 
all values of p. For low temperatures, on the other hand, 
the fit becomes worse (lower inset of Fig. 4). Similar 
results were presented in the experimental data of |l3| . 
The reason for this is that at low enough temperatures 
tunneling into the SCIs is suppressed, except very close 
to the SIT. Since the resistance through normal areas is 
activated not by E c but rather by W, the slope of the Ar- 
rhenius plot changes, as is also evident in the experiment. 
We thus predict that the low temperature behavior near 
B c will be the same as in higher temperatures in strong 
magnetic fields B » B max . 
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FIG. 4: The activation energy To as extracted from fitting the 
resistance to an activation behavior, obtained from the numerical 
calculation and from the experimental data of \\ A (inset). Lower 
inset : an Arrhenius plot of the resistance as a function of temper- 
ature for p = p m ax = 0.5. A deviation from an activated behavior 
is clearly seen. 

The fact that the temperature dependence is not a pure 
activation in the whole magnetic field and temperature 
range is crucial to the observation of such a peak in the 
MR. If one assumes solely a magnetic field dependent 
activation energy, Tq(B), then for different temperatures 
Ti and T2 the ratio 

log(R(B,Tl)/Ro) = Tx/Tp^) = J\ 
log{R{B,T2)/R ) T 2 /T (B) T 2 ' U 

would bc independent of field, in contrast with the ex- 
perimental observation. 

The amount of disorder affects our model in several 
ways. First, the width of the energy distribution W in- 
creases, though this has a minor effect on the behavior 
near the MR peak. Second the initial concentration of 
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normal islands, po, increases with disorder. If one as- 
sumes p(B) ~ po + aB x (where a is some constant and 
x is probably equal to 2 for small fields, as the system 
is symmetric under reversal of the magnetic field direc- 
tion), then B max = (p max -p Q ) 1/x /a, namely B max de- 
creases with increasing disorder. Lastly, the typical size 
of the SCIs decreases, leading to an enhancement of the 
Coulomb charging energy and an exponential increase of 
J? max . These latter two points are consistent with the 
experimental observations |T3 |. though quantitative pre- 
dictions require a detailed study of the distribution of the 
SCIs with disorder and magnetic field (currently under- 
way). 

As mentioned in 14], at very strong fields (up to 
30 T) the resistance of some samples saturates at val- 
ues somewhat larger than the resistance just above the 
(temperature-driven) superconducting transition. Our 
model provides a natural explanation for this. At strong 
magnetic fields all vestiges of SC correlations are gone, 
and one is left with an insulator (which is hardly affected 
by magnetic field) at a lower temperature, and thus with 
higher resistance. 

The theory presented here may also account for the 
MR in the presence of a tilted magnetic field 19]. Tak- 
ing into account the finite width of the sample, a tilted 
magnetic field suppresses the SCIs both in-plane and in 
the "thin" direction. Due to the finite width, the elec- 
tron trajectories are now 3-dimensional, and thus changes 
in the SCI size due to parallel magnetic field affects the 
resistance in a similar way as described above, leading 
to the same non-monotonic MR for parallel fields. This 
also explains the shift in the location of MR peak towards 
higher fields for parallel fields ■ 

To summarize, we have demonstrated that competition 
between normal electron and Cooper-pair transport, gen- 
erated by Coulomb blockade of superconducting islands 
and driven by perpendicular magnetic field, may yield 
non-monotonic magneto-resistance. This is accompanied 
by a change in the temperature dependence of the resis- 
tance, resulting from a crossover from nearest-neighbor 
hopping in normal areas to tunneling into the SC islands. 
This crossover may lead to a change of several orders of 
magnitude in the resistance, as seen experimentally. 
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